Recombination of H3 Ions in the Afterglow of a He-Ar-H2 Plasma 
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Recombination of H|J" with electrons was studied in a low temperature plasma in helium. The 
plasma recombination rate is driven by two body, HJ+ e~, and three-body, e~ + He, processes 

with the rate coefficients 7.5 x 10~* cm^s"-^ and 2.8 x 10~^^ cm'^s"-^ correspondingly at 260 K. The 
two-body rate coefficient is in excellent agreement with results from storage ring experiments and 
theoretical calculations. We suggest that the three-body recombination involves formation of highly 
excited Rydberg neutral H3 followed by an I- or m- changing collision with He. Plasma electron 
spectroscopy indicates the presence of H3. 

PACS numbers: 34.80.Lx, 34.80.Kw, 52.72.+V 
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Recombination of with electrons has been the sub- 
ject of many experimental and theoretical studies for 
nearly 50 years. Because of the high hydrogen abundance 
in the universe, Hj^ plays a key role in interstellar molec- 
ular clouds and planetary atmospheres It is also an 
important ion in many laboratory hydrogen-containing 
plasmas and in technological applications, including cold 
regions of fusion devices. Despite enormous effort, the 
process of recombination is still controversial (see the 
critical discussion in Ref. 0). 

Early theoretical studies suggested that the re- 
combination rate is slow, although experiment gave a 
much higher recombination rate. Inclusion of non- 
Born-Oppenhcimer Jahn- Teller coupling significantly im- 
proved the theory and [3, 4] gave a high rate of the disso- 
ciative recombination (DR) with the rate coefficient (RC) 
a = 7.8 X 10~^ cm'^s"^ at 300 K. Recent storage ring ex- 
periments d, [1] with cold ion sources obtained a DR rate 
for 300 K in good agreement with theory. Discrepancies 
remain just in connection with para and ortho H;^ DR. 
The disagreement with plasma experiments has not been 
solved until now 0, [1| . 

Essentially two different types of the experiments have 
measured the DR RCs: Colliding beam experiments [Q], 
where binary cross sections are measured, and afterglow 
experiments where one monitors the decay of the H;^- 
dominated plasma and the recombination RC a{T) ver- 
sus temperature is obtained. The essential difference is 
that in the beam experiments only binary collisions are 
important, whereas in a plasma many-body processes are 
possible and can play an important role. 

In this article we present new experimental data for 
recombination in an H;^-dominated plasma. The new 
results and the analysis reveal the importance of three- 
body processes in plasma recombination and finally rec- 
oncile the theoretical and experimental storage ring re- 
sults with the plasma experiments. 

The first experiment discussed here is the advance in- 
tegrated stationary afterglow experiment (AISA), where 



a plasma was generated by pulsed microwave discharge 
in a IIe/Ar/Il2 mixture f^,^. The presence of He/Ar/H2 
is essential because H^t is formed in the sequence of ion- 
molecule reactions [3l|. After the microwaves shut off, 
the decay of the cold dominated plasma is monitored 
over 60 ms with a Langmuir probe. From this plasma de- 
cays, recombination RCs as small as 5 x 10~^ cm^^s"^ can 
be measured with AISA. In AISA we studied the depen- 
dence of recombination RC on the partial pressure of H2. 
Examples of the data obtained are plotted in Fig. [T] The 
plasma recombination RC acff is extracted under the as- 
sumptions that the plasma is quasi-neutral and that the 
decay is dominated by processes described by the balance 
equation: 



-acSUen^ 



(1) 



where rie and n+ are the electron and ion densities, re- 
spectively. The AISA data plotted in Fig. [T] show that 
at low hydrogen densities, acff increases with [H2] until 
it reaches "saturation" at [H2] ~ 2 x lO^^cni^^. 

To verify that in our plasma experiments ions other 
than H;^ do not play an important role in recombination 
(see discussion in [l3|) we built another stationary after- 
glow experiment - Test Discharge Tube equipped with 
a near-infrared Cavity Ring Down Spectrometer, TDT- 
CRDS 11]. In this experiment the decay of H^(u — 0) 
was monitored during the afterglow (after microwaves 
shut off) and Ofeff was obtained. From the Doppler broad- 
ening of the absorption line the kinetic temperature of 
ions was determined [12i] . High electron (and ion) den- 
sity was used in the TDT-CRDS experiment in order to 
have sufficient signal-to-noise ratio. The results of TDT- 
CRDS obtained in the "saturated" region agree well with 
the AISA data. 

To obtain even higher control over the plasma param- 
eters in a recombination-dominated afterglow we have 
adapted a flowing afterglow apparatus (FALP variant 
[1^, see Fig. [5]) to study the slow recombination pro- 
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FIG. 1: The variation of recombination RC with H2 density 
for two different He pressures. The data were obtained in 
AISA and FALP experiments. The theoretical value for 
the binary RC is indicated_|3|, |4| , as is the current value from 
storage-ring experiments 



In this FALP, the 
plasma is formed upstream in the microwave discharge. 
To remove metastable helium by Penning ionization, Ar 
is added through port Pi. In the He/Ar mixture an Ar+ 
dominated plasma is formed 0, and carried along 
the flow tube for 30 ms. When the plasma is already 
cold, H2 is introduced via port P2, which forms an H;^ 
dominated plasma 0] . The decay of the dominated 
plasma is monitored by Langmuir probe along the flow 
tube for another 60 ms. The coefficients aes obtained by 
FALP for different [H2] at 1600 Pa and 250 K are plotted 
in Fig. [T] Again, a dependence of acff on [H2] was ob- 
served. It is also clear that ace depends on He pressure, 
suggesting that the recombination in a plasma is not just 
a binary process. Note the small increase of aos at [H2] 
approaching 2 x lO^^cm"^, this is due to the formation 
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ions in three-body association reaction 
He — > H^ -|- He, and their immediate recombination 
with electrons. This process was studied separately us- 
ing FALP and TDT-CRDS and we obtained very clear 
results an d ag reement with thermodynamic data on li^ 
formation On the base of these studies we can con- 
clude, that influence of H^ formation on measured actf 
can be excluded at 260 K for [H2] < 
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2x10 cm ^ < [H2] < 1x10 cm ^. In these conditions, 
ttoff depends only on the He pressure. Fig. [3] plots the 
values of aofr obtained in AISA, FALP and TDT-CRDS 
experiments in the "saturated" region at 260 K for vari- 
ous He densities. Fig. [3] also includes the theoretical RC 
[l9| and results obtained in other laboratories where a 
He/H2 mixture was used; in cases where their data were 
obtained at different temperature (but close to 260 K) 
the corresponding value aeff(260 K) was calculated as- 
suming a T~^'^ dependence. The following observations 
are relevant to interpreting the plotted data: 
(i) Pittsburgh group, stationary afterglow [l3l. The mea- 
surement was in the saturation region at 18.7 Torr and 
300 K. Qfoff = 2.3x10"^ cm^s^i (Fig. 2ofRef. (H). (ii) 
Smith and Spanel (Fig. 4 of Ref. [lOj). At the beginning 
of the decay aoff (300 K) ^ SxlQ-^cm^s"^ at 2 Torr. For 
sufficiently long decay time the recombination coefficient 
becomes very small, (iii) Laube et ai, FALP at 0.5 Torr. 
aeff(300 K) = 7.8 x lO^^cm^s^^ (abstract and Fig. 3 of 
Ref. 0). (iv) Adams and Smith, FALP. Fig 2. of Ref. 
[l6t shows the decay curve with a fast decay at the begin- 
ning and a slower decay towards the end of the flow tube. 
Interpretation for this behavior will be discussed else- 
where. Here, we use the rate constant of the fast decay, 
aoff(300 K) = 7 X lO^^^cm^s"! at 1.2 Torr. (v) Storage 
ring with cold ion sources @, il- (vi) Our lab^DT- 
CRDS, present data, and data from Refs. [ll|, El, [13 ■ 
These are the only measurements made in He/Ar/H2 
mixture using absorption spectroscopy. The decay was 
measured in relatively early afterglow, at 50 — 300 ^s and 
electron densities Ueo ~ 5x lO^^cm"'^. (vii) Our lab - sta- 
tionary afterglow jAISA and VT-AISA, present data and 
data from Refs. 0,[|. (viii) Our lab - FALP (two dif- 
ferent constructions), present study and data from Refs. 
0, [3 ■ (ix) Theoretical value [HI . Experimental data 
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FIG. 2: Principle of the FALP method. Buffer gas (He) flows 
from the discharge region towards a pump (right side). Re- 
actants are added via ports Pi and P2. The relation between 
time (t) and position is given by buffer gas velocity. 
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FIG. 3: The dependence of effective recombination RC Qeff 
on He number density. A straight line is plotted through the 
AISA and Pittsburgh data points. 



Region of high H2 density. In order to under- 
stand the mechanism of recombination we concentrate 
first on the "saturated" region, i.e. on data obtained with 



from 10 very different plasma experiments measured over 
24 years suggest that effective acs depends linearly on the 
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He density 

acff = + i^Hc[He] (2) 

with the coefficient of three-body recombination ~ 
2.8 X 10~^^cm^s~^. An extrapolation for [He]— >0 gives 
the value a°ff(260K) = 7.5 x lO^^cm^s"^ This value 
agrees well with the Ho" DR RC obtained in several stor- 
age ring experiments 

and theory 

The observed linear dependence of OcS as a function 
of the helium density can be attributed to an additional 
channel of recombination of 11^ , which is effective only 
in presence of He and can be called He-assisted recombi- 
nation: 

r H2 + H + He or 
H+ + + He < H + H + H + He or (3) 
I H^+Hc 

with the three-body RC i^Tne- The first two channels 
in Eq. [3] can be called He-assisted DR. The molecules 
H3 have a finite lifetime and can be relatively stable only 
in an excited electronic state. Such metastable molecules 
play am important role at low densities of background H2 
gas (see below). Previously, three-body processes have 
been observed and described by Thomson and later Bates 
[20| . Typical values of three-body RCs for He as a third 
body are 10~^^cm®s"^ (see refs. [13, Ull and references 
therein) . Processes with such RCs can contribute only at 
pressures higher than 1 atm. The process we observed is 
more efficient by factor of 100. It is however possible that 
the three-body process of Eq. ^ proceeds through an 
intermediate state of Hg formed during binary Hg -|- e~ 
collisions. If there is no external perturber (such as He 
or H2), the lifetime of such Hg is mainly determined by 
autoionization: the autoionization probability is much 
larger than the predissociation. Typical lifetimes of elec- 
tronic p-states of Hg (i.e., those most likely to be formed 
in Hg" -I- e~ collisions) are = 10 — 100 ps. If, dur- 
ing this time, Hg collides with a He atom followed by a 
change of the electronic state of Hg, (say, Z-changing Hg 
+ He collision) the new state should have a much longer 
lifetime with respect to autoionization. These states may 
be stable or not with respect to predissociation. We use 
the notation Hg (in contrast to Hg) in Eq. [3] and below to 
label the states stable with respect to predissociation and 
autoionization. It is possible to estimate the three-body 
RC k* for formation of such states. 

The number N* of Hg molecules created in 1 cm'^ per 1 
second is fc; [He] [Hg] , where ki is the RC for the /-changing 
collisions. Assuming that Hg are constantly being created 
in H;^ + e^ collisions and destroyed by autoionization, 
their density [Hg] can be estimated as [Hg] = [H;^]P, 
where P is the fraction of time that an ion H^ is being in 
the "form" of Hg. Thus, P = Ta/Tc, where Tc — l/(a*ne) 
is the time between HJ -|- e^ collisions, a* is the RC for 
H* formation. Therefore, TV* = fci[He][H+]TQ a* rig, and 
the three-body RC is estimated as k* = kiTaa* . 

The RC a* is larger than the DR RC a^R = 7.5 x 
lO-^cm^s-i (T = 260 K) Q by about a factor 10, 



i.e. a* ^ 8 X lO-'^cm^s'^ The RC k for the /- 
changing collisions can be estimated from experimental 
and theoretical data for Na*+He /-changing collisions 
[13]: ki ~ 2.3 X lO-^cm^s-^ If we take Ta ^ 20 ps, 
the estimated value for fc* is 4 x 10~^^ cm^s~^. 

Region of low H2 density. When the density of 
the H2 background gas is smaller than 2 x 10^^ cm~^ the 
AISA and FALP experiments show almost the linear de- 
pendence of effective aes on [H2]. This can be explained 
by the presence of long-living states of the neutral Hg 
molecule that are formed in the decaying plasma (the 
third channel in Eq. [3]) and can only be destroyed in 
collisions with H2 molecules. Such metastable molecules 
Hg can be formed in the /-changing collisions Hg 4- He 
discussed above. Because they can only be destroyed by 
H2, at low H2 densities such Hg molecules play a role 
of the "storage reservoir" for the H^ ions. Slow release 
of the H;^ ions from the Hg molecules makes measured 
effective acff to be dependent on the H2 density. 

Our calculation shows that the molecules Hg (having 
relatively short lifetime 10— 100 ps) are formed mostly in 
highly excited Rydberg np states with n '--^ 20 — 70. This 
is because the second ionization threshold for para-H;^ 
is 170 cm""^ above the lowest ionization threshold. This 
difference happens to coincide approximately with the 
temperature of the experiment T — 260 K. At energies 
below the second ionization threshold, the Hg system has 
an infinite number of p-state resonances with widths com- 
parable to the splitting between resonances. The proba- 
bility of capture to one of these resonances is comparable 
to unity ^]. Thus, if Hg is formed from Hg in /-changing 
collisions with He, they are formed in highly excited Ry- 
dberg states that should have very long lifetimes since 
the lifetime scales as n^. In fact, at even lower electron 
temperatures ^ 100 K, the binary DR RC is predicted to 
be significantly smaller [l9j for ortho-H;^ than for para- 
H;^. This should partially compensate for the "storage 
effect" at lower temperatures. The present experimental 
data suggest that such metastable molecules are indeed 
present in the decaying plasma, as is shown below. 

Thus, the presence of He and H2 influences the over- 
all process of plasma recombination: When collides 
with an electron, a helium atom can speed up the disso- 
ciation of Hg or help to create metastable molecules Hg . 
If there are many H2 diatomic molecules such metastable 
molecules are destroyed. It is the H2- "saturated" re- 
gion. In this region, the recombination RC is described 
by Eq. ([2]). When the density of H2 is smafl, the Hg 
molecules are accumulated and play a role of a tempo- 
rary storage for H;^ ions. For this region the effective 
recombination rate agff depends on [H2]. 

If Hg exist and if collisions with H2 influence the pro- 
cess of recombination, then if another reactant will be 
added with H2 we should observe the reaction. There 
are not many possibilities; if Hg has to be also produced 
Xe is probably the only choice. After addition of Xe a 
new channel is open: Hg -I- Xe — > XeH"*" -I- e^^^ -I- H2 . If 
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FIG. 4: EEDF measured in PALP experiment in 
He/Ar/H2/Xe decaying plasma: EEDFs in several times of 
the plasma decay are depicted. The data obtained have been 
fitted to a Maxwellian distribution (full lines) and the EEDF 
are normalized. For comparison we also show the EEDF (open 
circles) measured in an Ar"^ dominated plasma without addi- 
tion of the H2/Xe mixture. 

H3 is highly excited, with energy close to ionization, then 
this reaction is exothermic and the ejected electron can 
have energy up to 0.75 eV. Fast electrons with energy 
0.75 eV produced in cold plasma with electron tempera- 
ture Te ~ 0.03 eV can be in principle detected by mea- 
suring the electron energy distribution function (EEDF) 
[23I . |2J|. This method is named plasma electron spec- 
troscopy (PES) [25i, 261 ■ the experiment the discharge 
was in He, then Ar was added via port Pi, and 30 ms 
later, Xe and H2 were added through port P2 considered 
as i = 0. The densities of H2 and Xc were made simi- 
lar, in order to have comparable probabilities of reaction 



with Hg. 

Measured EEDFs in He/Ar/H2/Xe afterglow plasma 
are depicted in Fig. 01 The measured EEDFs clearly 
indicate a production of fast electrons with energy s ~ 
0.75 eV in the plasma, as expected for the Xe reaction 
process. The method is not yet sufficiently elaborated to 
account for potential superelastic collisions of electrons 
with internally excited ions (e.g. XeH"*") and to give 
quantitative results. Qualitatively, the measured EEDF 
evolution indicates that at low decay time, when the re- 
combination is fast, the production of energetic electrons 
is also fast; at a longer decay time the production be- 
comes slower and EEDF approaches a Maxwellian shape. 

Analysis of the present and several previous experi- 
ments involving an H;^-dominated plasma provides firm 
evidence that the rate of plasma recombination strongly 
depends on the density of the background helium. More 
importantly, the analysis suggests that at the limit of 
zero He density, the thermal RC for binary recombina- 
tion in all considered plasma experiments is in agree- 
ment with the recent storage ring experiments and with 
theory. Thus, the 50-year old problem of disagreement 
for the binary DR RC at low energies between differ- 
ent experiments and theory appears to be resolved. Fi- 
nally, the present experimental data suggests that long- 
lived H3 molecules exist in the Hj^-dominated plasma 
and play an important role in the plasma recombination. 
Acknowledgements: This work is a part of the research 
plan MSM 0021620834 financed by the Ministry of Edu- 
cation of the Czech Republic and partly was supported 
by GACR (202/05/P095, 205/05/0390, 202/03/H162, 
202/07/0495) by GAUK 53607 and GAUK 124707 and 
by the National Science Foundation, Grants No. PHY- 
0427460, PHY-0427376. 



[1] T. Oka, Phil. Trans. R. Soc. Lond. A 358, 2363 (2000). 

[2] R. Johnsen, J. Phys.: Conf Series 4, 83 (2005). 

[3] V. Kokoouline, C. H. Greene, B. D. Esry, Nature (Lon- 
don) 412, 891 (2001). 

[4] V. Kokoouline and C. H. Greene, Phys. Rev. A 68, 
012703 (2003). 

[5] B. J. McCaU et al, Nature 422, 500 (2003); B. J. McCaU 

et al. Phys. Rev. A 70, 052716 (2004). 
[6] H. Kreckel et al, Phys. Rev. Lett. 95, 263201 (2005). 
[7] R. Plasil et al, Int. J. Mass Spectrom. 218, 105 (2002). 
[8] J. Glosik et al, Chem. Phys. Lett. 331, 209 (2000); 

J. Glosik et al. J. Phys.B 34, L485 (2001). 
[9] A. Wolf et al, Phil. Trans. R. Soc. A 364, 2981 (2006). 
[10] D. Smith and P. Spanel, Int. J. of Mass Spectr. Ion Proc. 

129, 163 (1993). 
[11] P. Macko et al, Int. J. Mass Spectr. 233, 299 (2004). 
[12] P. Macko et al, Acta Physica Slovaca 54, 263 (2004). 
[13] J. Glosik et al, Plasma Source Sci. Technol. 12, S117 
(2003). 

[14] M. T. Leu, M. A. Biondy, and R. Johnsen, Phys. Rev. A 



8, 413 (1973). 
[15] S. Laube et al, J. Phys. B 31, 2111 (1998). 
[16] N. G. Adams, D. Smith, and E. Alge, J. Chem. Phys. 81, 

1778 (1984). 

[17] R. Plasil et al, J. Phys.: Conf. Series 4, 118 (2005). 
[18] A. Pysanenko et al, Czech. J. Phys. 52, D681 (2002). 
[19] S. Fonseca dos Santos, V. Kokoouline, and C. H. Greene, 

J. Chem. Phys. 127, 124309 (2007). 
[20] D. R. Bates and S. P. Khare, Proc. Phys. Soc. 85, 231 

(1964). 

[21] Y. S. Cao and R. Johnsen, J. Chem. Phys. 94, 5443 
(1991). 

[22] A. P. Hickman, Phys. Rev. A 18, 1339 (1978). 

[23] J. Glosik et al, Int. J. Mass Spectrom. 189, 103 (1999). 

[24] R. R. Arslanbekov, N. A. Khromov, and A. A. Kudrjavt- 
sev. Plasma Source Sci. Technol. 3, 528 (1994) 

[25] N. B. Kolokolov, A. A. Kudrjavtsev, and A. B. Blagoev, 
Physica Scripta 50, 371 (1994). 

[26] I. Korolov et al, Contrib. Plasma Phys., sent for publi- 
cation. 



